Creation of Kink and Antikink Pairs Forced By Radiation 
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The interaction between kink and radiation in nonlinear one-dimensional real scalar field is in- 
vestigated. The process of discrete vibrational mode excitation in (f> 4 model is considered. The role 
of this oscillations in creation of kink and antikink is discussed. Numerical results are presented as 
well as some attempts of analytical explanations. An intriguing fractal structure in parameter space 
dividing regions with creation and without is also presented. 
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INTRODUCTION 

Topological defects are usually compact, static solu- 
tions with finite energy of partial differential equations 
which properties are defined by field values at infinity. 
The simplest examples of topological defects are one- 
dimensional kinks ((p 4 or sine-Gordon equation). Less 
trivial examples (vortices, domain walls, monopoles etcj 
manifest themselves in various branches of physics Q , Q 
like particlephysics Q, condensed matter physics, cos- 
mology 0, Q and much more. 

Because of their stability and localized energy density 
some topological defects have similar properties as parti- 
cles (e.g. kinks in 1+1 dimensions or monopoles in 3+1 
dimensions). They can be created or annihilated as or- 
dinary particles. They can also interact with each other, 
with radiation or external force. |(| brings a very nice ex- 
ample of how topological defects accelerate in a constant 
external force. 

However, topological defects reveal some properties 
which make them very different from particles. In our 
previous paper Q we have presented an unexpected be- 
havior when the kink in <^> 4 was exposed to monochro- 
matic radiation of scalar field coming from one direction. 
The nonlinear character of this interaction resulted in 
"negative radiation pressure" , the kink, for most of the 
frequencies of incoming wave, instead of being pushed by 
radiation, was in fact being pulled toward the source of 
this radiation. Similar feature reveals the second most 
often used example, the sine-Gordon soliton. 

Another peculiarity of topological defects occurred in 
a two kink collision process in </> 4 theory || during which 
the kinks could be scattered back or annihilated. Let us 
stress out that similar process in sine-Gordon equation, 
because of its intcgrability is less interesting. If the SG 
solitons have enough energy they go through each other 
with nothing but a phase shift. When the energy is be- 
low some threshold the solution describe a bound state 
- a breather which is a perfectly periodic solution. (f) A 
has much more interesting structure of solutions. During 
collision the kinks loose some of their initial energy on 
behalf of radiation and when the loss is small the kinks 
reflect from each other. When the loss is large they glue 
to each other forming a bound state (so called oscillon 



|9()- The bound state also radiates and finally vanishes. 
An extraordinary thing about this process is that there 
exists no real threshold. The regions in phase space for 
one or the other scenario mix and form a fractal similar 
to a Cantor set. 

This is not an isolated example, brings a differ- 
ent example when a sine-Gordon soliton interacts with a 
oscillating impurity mode. They also observed a fractal 
structure. 

We found yet another example of fractal structure but 
in a quite different process. In (j> A the soliton possesses an 
internal degree of freedom - an oscillational mode (in jjj 
authors claim it is responsible for a peculiar behavior of 
kinks during collision) . When that mode is excited it os- 
cillates with a certain frequency. In linear approximation 
the mode oscillates infinitely long with constant ampli- 
tude. When one includes nonlinear correction one finds 
the mode couples to scattering spectrum and radiates re- 
ducing its amplitude and finally vanishes. But when the 
initial amplitude is large enough so that its energy is just 
a little above a mass of two kinks, a kink and antikink 
are produced. 

In present work we investigated an opposite process. 
We light radiation far away from the kink from both di- 
rections so that the kink would stand still. The radiation 
excites the oscillational mode and if the radiation am- 
plitude and frequency is suitable the oscillational mode 
breaks up to kink and antykink. Of course if radiation is 
large enough the kinks can be created even from vacuum, 
but because the oscillational mode gathers the energy 
our process seems to be much more efficient. One sur- 
prising result we have found is that in a plane radiation 
amplitude vs radiation frequency the border separating 
solutions with and without creation is also a fractal. We 
claim that the nonlinear coupling between the vibrational 
mode and radiation is responsible for the structure. 

Our paper is organized as follows. In the following 
section we give an introduction to the model discussed. 
We remind a spectral structure of linearization around 
the kink. The next section is devoted to excitation of the 
oscillations mode. We present our numerical results and 
attempts of theoretical explanation. We present a simple 
model reproducing qualitative results of the full model. 
Although on a precise level there are some discrepancies, 
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the general structure of solutions remains similar. The 
last section is as usual conclusions and discussion. 

We have encountered also a very interesting question 
of how one can simplify and reduce a system from infinite 
number of degrees of freedom to only few. 



THE MODEL 

Let us consider one-dimensional real scalar field obey- 
ing the equation in the form: 



4> - 4>" + U'{4>) = o, 



(i) 



where U(4>) is a potential with at least two equal minima 
(vacua) <p v (without loss of generality we can assume that 
U(4> v ) = 0) so that in these theories there exist static 
soliton solutions 4>s described by the equation 



4> a (x) : x - x Q = ± 



(2) 



In this section we will discuss two the most known models 
i.e. 4 (U(<f>) = \{(f) - l) 2 ) and sine-Gordon (U(<j>) = 
1 — cos((j))) equations. In these theories the static soliton 
solutions have the form: 



} (x) = ± tanh(x — x Q ) for </> 4 
,(x) = ±arctan(e-^- a: °)) for S - G 



(3) 



Because of the translational invariance we can substi- 
tute xq = 0. Let us add some small perturbation to the 
static kink solution ((f) = <f> s + £). If the potential can 
be expanded in a Taylor' series (one can do so for many 
systems with an exception of certain compactons |l2j). 

u'(4>) = u'fa + = u'fa) + c/"(0 s )e + N(t, 4>s) (4) 

where -/V(£, <f) s ) is a nonlinear part in £. We can write the 
equation for £: 

| - e" + V{x)S + N(£, x) = Z-L£ + N(£, x) = 0. (5) 

where 
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V(x) = U"(Mx)) 



cosh 2 x 



for S - G 



(i 



(6) 



cosh 2 x 



for 



We substitute £{t,x) = e luJkt r]k{x) and find the eigenval- 

-A + 



,2 

ues and the cigcnfunctions of the operator L - 



V(x). We can divide the spectra into three groups 

1. translational zero modes i] t = <f>' s : 4> s (x + 5x) 
4> s (x) + 6x<j)' s (x), Lu t = 0, 



2. discrete oscillational modes (there is none for S-G 
and only one for </> but in general there can be any 
finite number) 



rf*\x) 



sinh x 
cosh 2 x ' 



v/3, 



3. scattering modes: 



. e zkx (ik — tanhx) 
Vk(x) = <; eikx ^ 3tanh 2 x _ 1 _ p _ tanh cc) for 



for S - G 
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It is very significant that scattering modes in these mod- 
els have no reflection part. This is not a general feature. 
In fact it is quite rare, but for example spectra for V(x) 
given in a form: 



V(x) 



N{N + 1) 
cosh 2 x 



are reflectionless for all integer N. There are exactly N 
bounded modes. One of them is a translational mode 
and the rest of them are oscillational. 

The existence of the oscillational mode in the (j> A theory 
is responsible for its nonintegrability. During for exam- 
ple collisions this mode is excited and takes some of the 
initial energy. Then the energy, due to the nonlinearities, 
is transfered to radiation modes. 

Numerical results show that the radiation coming from 
infinity and hitting the kink can cause the creation of 
the pair kink and antikink. This process is much eas- 
ier for <fi 4 . Even small amplitudes can force the process. 
In sine-Gordon's case the process occurs for amplitudes 
compared to the amplitudes when the creation can occur 
even far away form kink in pure radiation. The explana- 
tion is quite simple. In </> 4 there exist an oscillation mode 
which can be excited by the radiation. Manton in- 
vestigated the creation of kink and antikink caused only 
by a discrete mode. When the energy contained in this 
mode was a little above then the energy of two kinks a 
pair of kink and antikink was created and two kinks were 
radiated to infinity in opposite directions and antikink 
remained. When the energy was not large enough the 
discrete mode oscillated with a decreasing amplitude and 
the energy was radiated to infinity because of the cou- 
pling to the scattering modes. In the following section 
we will discuss an opposite process. Radiation coming 
from infinity will couple to the oscillating mode. 



EXCITATION OF THE OSCILLATING MODE 
AND THE CREATION PROCESS IN 4 MODEL 

Numerical results 

We used two different sets of initial and boundary con- 
ditions to determine the coupling between radiation and 
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FIG. 1: The excitation of the discrete mode for A = 0.2 and 
ui = 3.5. 
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FIG. 2: The excitation of the discrete mode for A — 0.3 and 
lj = 3.5. 



the oscillating mode. In both cases we considered only 
antisymmetric problems, we evaluated our system only 
on one half of the x-axis and posed the boundary condi- 
tion <fi(t, 0) = 0. We could do so because the symmetrical 
radiation does not couple with an oscillational mode and 
therefore gives no contribution. The symmetrical radia- 
tion couples however with a translational mode and the 
kink as a whole would start moving in one direction and 
would be a difficult object to study. 

In the first case we had initial conditions 0(0, x) = 
(j> s {x). The boundary conditions were (j)(0,x) = and 
4>(t,L) = ^As'mojt (where L was large in comparison 
to the size of the kink, usually L = 200). We use | 
because when the waves come from both sides of the 
kink, they add and produce a standing wave with twice 
as big amplitude. Because of the condition <fi(t, 0) = 
and the symmetry of the equation <f> — > —<fi, x — > —x we 
actually studied the system with boundary conditions 
4>(t,±L) = ±7}Asmujt. These conditions describe the 
kink in radiation coming from both sides from very far 
away. In order to measure the excitation of the discrete 
mode we calculated the value (a projection on a discrete 
mode): 

/oo 
dxr) d (x) [<f)(t,x) - <f> s (x)] 

AdW = — Too ■ (7) 

/ dxr)j(x) 

J — oo 

Because the radiation is orthogonal to r\d there is no 
contribution to Ad{t) from any other modes (at least in 
linear approximation). We measured Ad for different w 
and A. 

In figures Hl3l we have sketched three examples of Ad 
vs time for ui = 3.5 and for A equals 0.2, 0.3 and 0.5. 
For all Figures we can see that the reaction of the os- 
cillating mode is retarded by about 240. Because the 
radiation is "switched on" at x = L it needs time 
t = ||i = ^ ~ 243 to get to the kink. 
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FIG. 3: The excitation of the discrete mode for A — 0.5 and 
uj = 3.5. We can see the creation of kink-antikink pairs for 
t w 500, t w 560 and t w 590. 



If the amplitude is small (Fig. ^) Ad oscillates with eigen- 
frequency Ud- These oscillations are shifted and modu- 
lated by another frequency. As we will show in the fol- 
lowing section this modulating frequency is oj m = u> — 2ujd 
(Fig. UJ. In our example oj m = 0.0359 and the period is 
175 so we can see a very good agreement. 

If the amplitude is larger the modulations arc with 
higher frequency, and the oscillations decay. This decay 
can be explained because as Manton showed in the 
oscillating mode couples to the scattering modes and ra- 
diates its energy (see also 0, 0). 
The next Figure Q is the most interesting. First we can 
see some chaotic oscillations and then around t = 500 Ad 
suddenly goes to about —1.5 and then oscillates around 
that level, and then goes back to (around t = 560). If 
the amplitude of Ad is large enough a pair of kink and 
antikink can be created and radiated in both directions, 
but the place of a kink in the middle is taken by an an- 
tikink remains. When we substitute <f>(x,t) = —<p s (x) in 
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FIG. 4: Fourier transform of A<j for w = 3.5 and A = 0.2. 
One can see four separated peaks for u = y3 ~ 1.73, 3.5, 7 
and 3.5 - ^3 « 1.76. 
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FIG. 6: Dependence 2 jp" 8 ^ upon logn;, . We can expect that 
a fractal dimension would be more than 1.6. 
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FIG. 5: Minima of Ad vs frequency u> and amplitude A of 
radiation coming from L = 200 . Dark spots Ad < — ^ w 
— 1.57 represent creation process. 



we obtain that Ad = — ? ~ —1.57. Figure [S] presents 
the minima of -i4d(£) as a dependence upon frequency w 
and amplitude A of incoming radiation. When Ad is less 
than —fa creation occurred (dark spots). One can see 
the border between creation and small oscillations is very 
complicated, perhaps even fractal. The resolution is not 
good enough to prove this proposition but the Figure 
may justify our hypothesis. We present the dependence 

D(logn fe ) = (8) 
log n b 

where n& is a number of boxes and I is a number of boxes 
containing the boundary. This value in limit rib ~~ * 00 
would give a true fractal dimension. Because we have 
only 2500 boxes we cannot tell the true dimension but 
it seems to be much larger than 1 somewhere around 
1.6 — 1.7 and that would justify that the border is fractal. 
Of course proving numerically that something is a fractal 
is very difficult and in our case is almost impossible. 



Although the case described above has a nice physical 
interpretation, because we switch on the light in one par- 
ticular moment, it is very difficult to explain it in details 
using a simplified method presented in the next section. 
The reason is that the simplified equations for Ad are de- 
rived for monochromatic wave coming from both sides. 
When we consider the boundary conditions described in 
the beginning of the present section we find that the radi- 
ation has a form of a very wide wavelet. Despite nonlin- 
earities in the equation the wavelet is always a superpo- 
sition of many monochromatic waves. Although we can 
assume that after a long time we have almost a single 
frequency standing wave there is still a problem of initial 
conditions for Ad(to). We found that the evolution of 
Ad{t) highly depends on the initial conditions. We can- 
not tell whether Ad would oscillate or jump to — f (and 
force the creation of soliton pairs). 

It is much easier to consider the Cauchy problem for con- 
ditions: £(£ = 0, x) = Ark(t = 0, #), where rk(t,x) is a 
real combination of scattering modes rj±k in a form of a 
standing wave r/c(i, x) = hk(x) cos tot where: 

(3 tanh 2 x — 1 — k 2 ) sin kx — 3 tanh x cos kx 



(9) 

The denominator was introduced just for convenience so 
that the amplitude of the wave far away from the kink 
would be 1 . Basically the structure of the solutions look 
the same but for different parameters A and the oscil- 
lations begin for t = 0. We have plotted an analogical 
Figure to 03 (Figure EJ ■ We can see that the border be- 
tween the two types of solution (with a creation (dark 
spots) and without) also seems to be fractal. 

If we make a similar figure using the longer evolution 
the border would be shifted in the direction of smaller A 
for obvious reasons. We tried to evolve our system as long 
as the border does not change significantly with longer 
time. Unfortunately this procedure was very difficult to 
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FIG. 7: Minima of Ad vs frequency to and amplitude A of 
radiation in case of a standing wave. 



apply and we cannot be certain if the border presented 
in Figured is a true limit, but we think it is a quite good 
approximation nevertheless. 



Theory 

Let us consider a simplified field equation for £ in the 
</> 4 model: 

e+£e+6^e 2 = o. (io) 

We have neglected only the cubic term. 
Let £ has a form 

£(t,x) = Ar k (t,x) + A d (t)rj d (x) + r){t,x), (11) 

where 77 is orthogonal to r\d- It is consistent to assume 
that Ad and 77 are both of order 0(A 2 ). After substitu- 
tion i|ll|) into HI U|l we obtain (in the second order in A): 

(A d + ujAd) m + V + Lri+ 
+6<f> s (\A 2 h 2 (1 + cos 2u)t) + A d r\\ + r/ 2 + 

+2Adi] d Ahk cosw/j + 2Adrjrjd + 2 Ah^j] cos ujt 



where 



a (k) = 3^ 



O 

dx cj) s (x)h 2 k (x)r]d{x) 



dx r) 2 d (x) 



(13) 



(14) 



(3(k) = 12 



dx (j} s (x)hk(x)r] 2 (x) 



(15) 



7 = 6 



dx r] 2 d(x) 

J — OO 

/OO 
dx 4> s {x)rf d (x) 
-OO 



dx T]j(x) 



9?r 
16' 



(16) 



The first two integrals can also be calculated analyti- 
cally: 



a(k) 



9tt 
64N< 



(8k 4 + 34k 2 + 17) 1 



1 



cosh kir 



0(k) = 3irt 



k 4 + 2k 2 -8 
8ATsinh^ 



(17) 



(18) 



(12) 

The first correction to A d originated from 77 would be 
of order 0(A d ) and now we will neglect this term. Since 
rj is orthogonal to 77^ we can multiply both sides of this 
equation by r\d and integrate (similar procedure was ap- 
plied in 0). The projected equation yields: 

A d +uo 2 d A d +A 2 a{k) (1 + cos 2uot)+(3(k)AA d cos ut+-f A 2 d = 0, 



where N = yj(k 2 + A){k 2 + 1). 

Equation (|13f) is very similar to Mathieu's equation 
but with extra driving force and nonlinear term. This 
is just a simplification because the oscillating mode cou- 
ples to the scattering modes and is a source of radiation 
itself(|l3, 0). This radiation carries away the energy 
from the discrete mode causing the decay of the mode 
and that means we should include a damping term in 
()13JI . We could predict the loss of this energy in the 
similar way as we did in but in this case there are 
some problems. The nonlinear part in eq. I)12|l contains 
sources in the form of rkA d rj d which means that we do 
not know how much energy is taken from the oscillational 
mode and how much from original radiation. Finally we 
do not know how this loss would decrease amplitudes for 
different frequencies of the oscillational mode. This is 
why we did not add the damping term. 

Let us find a solution of the equation in a form (there 
are no free oscillations of the discrete mode 77^ i.e. there 
are no OA and OA 1 terms): 



AfA 2 



AfA z + AfA A + --- 



with a condition that A is very small. The first equation 
is of order 0{A 2 ): 



A d 2) + ujjA d 2) + a(l + cos2cji) = 0. 



(19) 



This equation is a simple forced harmonic oscillator equa- 
tion and we can solve it very easily: 



ALU 2 



cos 2u)t + B (Z) cos(uo d t + 5) . 

(20) 



We can also apply initial conditions such as Ad(0) = 
0, Ad(0) = and from that we can calculate and 
6. Let us assume (5 = 0. We can write the next order 
equation: 



.4 



(3) 



(3) 



t3A { p cos ujt 



0, 



(21) 



and solve it as well. The inhomogcncous term is equal 
to: 



(3A, 2) cos ujt = a(3[- 



cos cut cos 3u)t + cos tut 



+^B^j3(cos(uJd + u)t + cos(w - uj d )t). 

(22) 

This is a source for harmonic equation. The so- 
lution would oscillate with the following frequencies: 
ui, 3lo, LOd + u> — uid (compare with Figure |3J. Of course 
there is a special case when to = 2uid- We have then the 
source with resonance frequency and the linear growth 
of the amplitude and therefore the small amplitude ap- 
proximation fails. In fact this approximation also fails if 
iJ 1 — Aujj is very small, because in the solution there is 
a term 
small. 



u> 2 -iu< 



and although A is small Ad is no longer 



The most interesting question now is how big the am- 
plitude Ad can in be. We neglected some terms which 
are important to solve this problem. First of all this is 
only the solution in 0(A 3 ) order. Numerical solutions 
of the equation H3fl show that the most important role 
plays the quadratic term in Ad- We also neglected the 
cubic term and of course the field r\ which is the field 
of radiation escaping from this oscillating mode. We do 
not know how small A should be around the resonance 
frequency so that the oscillating mode amplitude would 
remain small during the whole evolution. 

Let us now focus on the equation for u> = 2tOd- Al- 
though we can find the precise values of the coefficients 
a(k) and f3{k) for k = 2v2 for this particular model we 
think it is better to examine the whole equation for any 
set of coefficients. For different models there should be 
different values but the structure of the equation remains 
the same. First we simplify the equation in order to get 
rid of too many parameters. We introduce 



(23) 



and our equation takes a form: 

w + w + w 2 + giwcos2r + g 2 (l + cos4r) = 0, (24) 

where ' denotes -4- Solutions of the above equation can be 
divided into three groups. In the first group the solution 
oscillates around some xq with frequencies w, 2lo, ■ ■ ■ and 
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FIG. 8: Solution of the eq. (ZQ for g 1 = 0.2, g 2 = 0.03. Low 
frequency modulation with large amplitude. 
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FIG. 9: Solution of the eq. (3J for gi = 0.1, p 2 = 0.03. High 
frequency modulation with small amplitude. 



is modulated with a very small frequency. The solutions 
have sometimes quite large amplitude, and the maximum 
is well above (Figure |SJ. 

The second group are solutions that oscillate between 
some minimum and 0. The modulation frequency can be 
both low and high (Figure |5J). 

The third group of these solutions are solutions which 
explode in a certain time r cr , they are singular (Figure 

us}. 

For large modulus of w we can neglect linear terms in 
w and focus on the equation: 



id + w 



0, 



(25) 



we can find the behavior of the solution for very large 
\w\: 



w 



-6(r - T a y 



Of course we cannot extrapolate these solution exactly to 
our partial equation for obvious reasons. First of all the 
amplitude in a field theory cannot grow unlimited, but 
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FIG. 10: Solution of the eq. 
Singularity around t = 57. 
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0.03. FIG. 12: Stable and unstable solutions on a plain (too, ii>o) for 
fixed pi = 0.5 and g2 = 0.12. 
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FIG. 11: Stable and unstable (dark) solutions on a plain 
(SiiS^) with initial conditions w(Q) = and w(0) = 0. 



FIG. 13: Minima of the solutions of the simplified equation 
in a plain of amplitude and frequency. 



in this case pairs of soliton and antisoliton are created. 
This process can also happen when the amplitude is large 
enough. 

There is also one more interesting observation accord- 
ing this simplified equation. The boundary between reg- 
ular and singular solutions on a plain g\ 1 g^ has a very 
complicated, presumably fractal, shape (Figure lTT|) . The 
same feature possesses the border between regular and 
singular solutions on a phase space of initial conditions 
with fixed 171 and gi (Figure ^J). 

In fact it is not surprising. When we take a closer 
look at the equation (|13fl we can see it possesses an in- 
variancc under discrete translation in time t — > t + T, 
where T = — . If we knew the map in phase space 
F: {A d (t),A d (t)) (A d (t + T),A d (t + T)) we could just 
iterate F n and find out whether for n — > 00 the solutions 
are singular or not. The same procedure one usually ap- 
plies when one wants to obtain a Julia or Mandelbrot sets 
some of which are one of the most known fractals. As 
a final result we present the Figure ITS! which shows the 



minima of the solutions of the simplified equation H3fl . 
Dark regions as usual represent the unstable solutions 
which certainly cross the level of — For those condi- 
tions we expect a pair of kink and antikink is produced. 

After comparing these results with the results obtained 
for full partial differential equation (Fig and EJ we can 
see that there are of course discrepancies but the the- 
oretical boundary lies not very far away from the true 
boundaries. It is actually easy to explain if we consider 
the energy which is radiated from the oscillational mode. 
In the theoretical picture we neglected this radiation. Of 
course we also did not include the cubic term and we 
neglected the orthogonal modes, but even in that over- 
simplified picture it is clear that something interesting 
happens. 

There is also one more interesting thing. Both in || 
and [lfij the simplified theories which were based upon 
collective coordinates reproduced the fractal structure 
which was wider than in field theory. In our case we used 
the projection onto the internal mode of the soliton and 
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our approach gave the structure which was much more 
narrow than in the partial differential equation. 

CONCLUSIONS AND DISCUSSION 

In the present paper we have presented the evolution 
of a kink with an external perturbation (radiation) in </> 4 
theory. We found that during the evolution an oscillating 
mode is being excited. We measured the excitation using 
numerical methods for two sets of initial and boundary 
conditions. Depending on the frequency of the radiation 
and its amplitude the excitation of the oscillational mode 
can remain small enough to be treated only as a defor- 
mation of the kink or can grow and finally can effect in a 
production of kink and antikink which are radiated into 
the spatial infinity. We tried to explain this process using 
a simplified method by projecting our field theory equa- 
tion onto the oscillational mode. T he obtained equation 
is a nonlinear ordinary equation which we cannot solve 
analytically, but we may study it numerically. This equa- 
tion possesses very rich space of solutions. Some of them 
are regular ones and some of them are singular. The 
border between those types of solutions is a very com- 
plicated both in a space of parameters and initial condi- 
tions. The solutions of the simplified equation are rela- 
tively good approximations of the solutions of the whole 
partial equation. The most interesting case of creation 
of kink-antikink pair reveals also a quite well agreement 
between field theory and our simplified theory. 
There are still many unanswered questions. First of them 
is why the border in the field theory is so complicated. 
We do not know how to add a damping term which we 
think is very important to solve the previous problem. 
Our simplified theory does not give the answers about 
the dynamics of the created defects. We do not know how 



much energy from the oscillational mode is converted into 
the kinetic energy of kinks, how much into radiation and 
how much energy remains in oscillating modes around all 
the defects. 
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